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Introduction 

, A quasi-optical macroscopic quantum phenomenon of birefringence, i.e., the effect of particle spin 
' rotation (oscillation) and spin dichroism accompanying the passage of high-energy particles with spin 
■ S > 1 through isotropic homogeneous matter was first described in [1, 2]. This effect is analogous to 
! the effect known in optics as the birefringence of light in an optically anisotropic medium associated 
with the dependence of the index of refraction on the state of light polarization, i.e., on the photon 
spin state. In contrast to light, whose wavelength is much longer than the distance between the 
atoms of matter, the de Broglie wavelength of a fast deuteron is much smaller than this distance. 
Nevertheless, according to [1, 2], in this case one can also introduce the spin-dependent index of 



Q^' refraction for particles. The deuteron birefringence effect occurs even in a homogeneous isotropic 
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medium and is due to the intrinsic anisotropy inherent in particles with spin S > 1. Moreover, 
CD I according to [1, 2], the magnitude of this effect even increases with growing particle energy, i.e., with 
diminishing wavelength. 

It was shown in [1, 2] that the spin dichroism described by the imaginary part of the index of 
refraction causes tensor polarization of the initially unpolarized deuteron beam. Recently, this effect 
00 ' was experimentally revealed in the energy regions 5 20 MeV [3, 4] and 3 GeV [5]. According to [4, 
' 5], the appearance of tensor polarization due to spin dichroism in particles with spin S > 1 can be 
used for designing the source of deuterons, which possess tensor polarization. According to [6, 7, 8], 
the birefringence effect should also be taken into account in experiments to search for deuteron EDM. 

In theoretical description of the birefringence phenomenon, it should be taken into account that the 
amplitude of forward scattering of a deuteron by a nucleus, which determines the index of refraction, 
includes the contribution due to Coulomb interaction alongside with that due to strong interaction. 
The Coulomb-nuclear interference correction to the rotation angle of the polarization vector and spin 
^ . dichroism of the deuteron was computed in [2, 9]. According to [9], the Coulomb-nuclear interference 
H I enables explanation of the change of sign of spin dichroism with changing deuteron energy in the 
range of 5 ^ 20 MeV. 

It should be noted that the interaction of deuterons and matter is accompanied not only by 
coherent scattering, which leads to the formation of a coherent wave, but also by incoherent scattering 
resulting in single and multiple scattering of particles in the target. In a kinetic equation for the density 
matrix, the characteristics of elastically scattered particles are described by the terms proportional 
to the squared absolute value of the total scattering amplitude (this amplitude is the sum of the 
Coulomb amplitude and the nuclear scattering amplitude modified by the Coulomb interaction) and 
as a result, in this case the terms describing the interference of the Coulomb and nuclear interactions 
should also be considered. In this connection, [10, 11] give the analysis of the correction due to the 
stated interference terms to the rotation of the polarization vector of protons passing through polarized 
matter for such target thicknesses and scattering angles at which the particle angular distribution is 
determined by the multiple Coulomb scattering and the single nuclear scattering in matter. It was 
shown in [11] that in the expression for the spin rotation angle of a high energy proton beam, the 
contribution of the Coulomb-nuclear interference to the amplitude of nuclear scattering at zero angle 
is compensated by the correction from the incoherent Coulomb-nuclear scattering. This statement 
was also formulated for scattering of particles with spin 1 [12]. 

In the present paper is shown that within the domain of applicability of the first Born approxima- 
tion, for the Coulomb amplitude the contributions to the rotation angle of the deuteron polarization 



vector coming from the Coulomb-nuclear interference are compensated only in the case when the 
detector provides a maximum coverage within a 47r solid angle geometry. If the detector registers the 
particles moving within a certain solid angle AO, <^ 47r with respect to the initial direction of the beam 
incidence, as it does in a real experiment, the interference terms stated above are not compensated. 
The magnitude of the contribution of the Coulomb-nuclear interference depends appreciably on the 
specific geometry of the experiment and is quite observable in the experiments on measuring the spin 
rotation angle of particles passing through the target. Multiple scattering of particles in the target 
leads to the fact that the magnitude of the contribution of the Coulomb-nuclear interference to the 
rotation angle of the beam polarization vector depends not only on the angle AO of the detector, but 
also on the mean square angle of multiple scattering. It is shown that if after the beam has passed 
the target, the magnitude of the mean square angle of multiple scattering is of the order of magnitude 
or greater than the squared diffraction angle of nuclear scattering, the Coulomb-nuclear interference 
can be neglected for any registered angle ACl of the detector. Otherwise, the contribution of the total 
Coulomb-nuclear interference to the rotation angle should be taken into account. 

1 Kinetic equation for the spin density matrix 

Let us consider the process of particle (beam of particles) passage through a target consisting of Nt 
particles interacting with one another by means of a certain potential U. 

Suppose that an incident particle has a rest mass m and spin Sd- The target consists of Nf bound 
particles with mass M and spin s. 

The Hamiltonian of the scattering system is written in the form: 

Nt 

Ht = Y,Ka + U, (1.1) 

a=l 

where Ka is the kinetic energy operator of particle a, U is the interaction energy of Nt particles of 
the scattcrer. 

The solution of the stationary Schrodinger equation 

Ht^j = W^^j (1.2) 

determines the possible values of the energy of the system and the corresponding set of wave 
functions. 

■^^ = -^^{Ri,si,...,Rn,sn), (1-3) 

where (Ra, Sq) are the spatial and spin coordinates of particle a in the scatterer. 

The operator of interaction between the incident particle and the scatterer is defined in terms of 

Nt 

V: V = ^ Va- Here Va describes the interaction between the beam particle and the target particle 

a=l 

a. The Hamiltonian of the whole system has a form: 

H = Kd + Ht + V, (1.4) 

where Ka is the kinetic energy operator of the beam particle d. 

To describe the process of particle d transmission through the target, find the density matrix 
p{t) of the system "incident particle -|- target". This density matrix satisfies the quantum Liouville 
equation: 

^^ = [H,pit)]■ (1-5) 
The solution of this equation can formally be written using the evolution operator U (t, to) 



p{t) = U{t,to)p{to)U+{t,to), 



(1.6) 



which is related to the exphcitly time-independent Hamiltonian of the system as U{t, to) = e~H^^*~*°^. 
The time moments to and t correspond to the state of the system before and after scattering of the 
incident particle by the a-th scattcrcr. 

Let us consider the target as a thermostat {Nf ^ 1). Then the statistical operator p of the 
system can be represented as a direct product p = Pd® Pt-, where pd is the density matrix of the 
incident particle, pt is the equilibrium density matrix of the medium. The spin density matrix p^ 
includes the elements diagonal and nondiagonal with respect to the momenta k of the incident particle: 
Pd = Pd{k, k'). However, the nondiagonal elements oscillate fast and after several collisions with target 
nuclei, in order to describe the process of multiple scattering, one can assume that the density matrix 
Pd{k,k') is diagonal [10, 13, 14, 15], i.e., Pd{k,k') = 6j^j^,pd{k). 

The time interval At during which the density matrix is diagonalized satisfies the inequality 
At 3> R/v, where R is the radius of action of the forces, v is the particle mean velocity in matter. 
This means that the evolution operator U{t,to) can be replaced by the Heisenberg's S-matrix = 

lim U{t, to), which relates the asymptotic states of the system before scattering to those after 

t— >-oo,to— >— oo 

scattering [10, 13, 14, 15]. The matrix elements of the S-matrix for a scattering system consisting of 
Nt particles are defined as follows: 



Nt 

Sba = ha - ^(27^)5(^6 - ^'^<-)ha . (1-7) 

a=l 

where Ea and Ei, are the total energies of the system before and after scattering, respectively; Ea = 
£k + W-y, El, = Ek' + W-y', £k and efe/ are the energies of the incident particle before and after the 
collision. In the momentum approximation, choose as Ta the scattering matrix of particle d interacting 
with a free particle a [16]. 

The limits to ^ — oo, t ^ oo are understood as the times when the incident particle is situated 
at distances larger than R, but the interval At = t — to is small as compared to l/v {I is the mean 
free path of a particle in matter). 

Denote the density matrix of the system at time t = to + At by p' , accordingly before scattering 
by the the a-th particle of the target, the density matrix will read p. Equation (1.6) can be rewritten 
for the density matrix of the incident particle, using the S-matrix defined in equation (1.7): 

p'a = Sp,5p5+, (1.8) 

where Sp^ means taking the trace over the states of the target. 

Write equality (1-8) for the diagonal momentum-space elements of the particle density matrix: 

At ^ - - Ai - - 

p'dik) = Pd{k) - i(27r)— Sp, Y,{k, -i\fa\k, ^)p{k, 7) + ^(27r)— Sp, ^ p{k, ^){k, ^\f+\k, 7)+ 

+(2,r)=— Sp, j;{t.7|t,|t',7'>/i(i',7')'5(eF -£» + WV - W.,){k' n'lf^lkn}- 

«i^=l k'j' 

In the momentum space, the matrix elements for the scattering matrix operator T have a form 
[17]: 

{k',P;,\fa\k,Pa) = {2'K)^5{k' + P'^ -k-Pa){k',P;,\Ta\k,Pa), (1-10) 

Pa denotes the momenta of the a-th scatterer, is the matrix of scattering by the momentum shell. 
Recall that Ta and Tq remain the operators with respect to spin variables. 

The states with a definite value of the momentum \k) are normalized according to equality {k'\k) = 

^ 6{k' — k), where V is the normalization volume. Thus, substitution of summation over all k by 



integration is made as follows y — )■ ^^2"^ Kronecker symbol and the Dirac function 5 are 

k 



related as J^^g, — > — Yj—6{k — k'). 
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Scattering processes that make the contribution to (1.9) can be both elastic and inelastic. Elastic 
scattering occurs in the absence of any excitation in the scatterer. Inelastic scattering accompanied 
by a small excitation in the scatterer is called "quasi-elastic" [16]. For such scattering the momentum 
q transferred to the a-th scatterer is g <C V 2MKa , q <C ^2MUa, q = k' — k; k is the momentum of 
the incident particle before the collision with the scatterer, k' is the momentum of the beam particle 
after the collision. 

When the momentum transferred to the scatterer a exceeds the momentum of the target particle 
in the initial bound state q 3> y/2MKa-, q 3> ^/2MTJ^, then the so-called "quasi- free" approximation 
can be used. In this approximation the excitation energy of the system in one collision is exactly equal 
to the recoil energy q^ j2M of a free particle of the target. Note that in summation over different 
nuclei of the target, sums of the form given below appear in the last term in (1.10) 



J d^Ri..J^RNe-''^('^--'^^^pt{Ri,...RN) = Ntpt+ 



a,/3=l 
Nt 



(1.11) 



+ Y j d^Ri..A^RNe-'^^^'^''-'^^^pt{Ri,...RN), 



where pt is the spin density matrix of the target nucleus. In deriving equation (1-11) it is assumed 
that the positions of target nuclei and their spin states are uncorrelated. 

Upon averaging over the spin states of nuclei in the target, the second term in (1.11) can be 
expressed in terms of the so-called particle pair distribution function [16], and vanishes when the 
transferred momentum q exceeds the magnitude inverse to the correlation radius r, i.e., q ^ [16]. 
For a non-crystalline target, the magnitude of the correlation radius r is of the order of the distance 
between nuclei. Consequently, the second term can only contribute to the kinetic equation at very 
small scattering angles -dgc ^ 1/kr, so it will be neglected in further consideration. 

The resulting expression for the density matrix can be obtained from (1.9) for ^'^ ^j^'^ ™ the form 
of the following integro-differential equation: 

dpd{k,t) 



-iVNtSTp^ (jik, k)p{k, t) - p{k, t)T^{k, kfj 



+ 



y3 r ) ' (1-12) 



^-^NtSp.J d'^k'S [eu-Sk' - f{kAk',-q)p{k',t)T^ik',-q;k,0). 



(2 

p{k) denotes the following dependnce: p{k) = pd{k;Sd) <8) Pt{St), where pt{St) is the spin density 
matrix of the target nucleus. 

Let us introduce the scattering amplitude F with matrix elements equal to [17]: 

F{k,k') = -MLv^T:(k,0;k',-q), (1.13) 

mM . , , , 

where Mr = — is the reduced mass. 

m + M 

Then equation (1.12) can be transformed into a form: 

^^4r^ = ^^Sp, (F{k,k)p{Kt)-p{UF^{k,k)) + 

r . . (1-14) 

^""'"^ J ^^ r,,./^- Fik,k')pik',t)F^ik',k), 



where N denotes the number of particles in the target per unit volume, ft' is the unit vector in the 
direction of the momentum k'. The absolute value of vector k' is determined from the equation: 

eu = ek' + - —. (1-1^) 

Note that when the condition m > M for the masses of the incident particles is fulfilled, the denomi- 
nator in the integrand of equation (1.14) vanishes for the value of the scattering angle of the incident 

particle equal to cos 9 = , with the absolute value of vector k' being equal to k\ -. 

m \ m + M 

Equation (1.14) simplifies when a particle (proton, deuteron, antiproton) passes through a target 
with nuclei whose mass is much larger than the mass of the incoming particle. In this case we can 
neglect the effect of the energy loss of the incident particle through scattering. So we can neglect the 
recoil energy (f /2M in the 5-function (1.12), (1.15). As a result, one obtains a simple kinetic equation 
describing the time and spin evolution of the incident particle as it passes through the target [10] : 

^MM. = ^sp, \F{k, k)pd{k, z) - pS, z)F+{k, k)] + iV-Sp, / d^-^,F{k, k'Ypik', t)F+{k', k), 

(1.16) 

where \k\ = \k'\. 

The first term on the right-hand side of (1.16), which describes refraction of particle in the target, 
can be represented as follows: 

F{0)p{k,t)-p{k,t)F+{0)= [i(#(0) + F+(0)),p(^,t)] +{i(F(0)-F+(0)),p(fc,i)}, (1.17) 

where [, ] is the commutator, { , } is the anticommutator. 

The part proportional to the commutator leads to the rotation of the polarization vector due to 
elastic coherent scattering (as a result of the refraction effect [10]); the anticommutator describes the 
reduction in the intensity and depolarization of the beam which has passed through the target. The 
last term in (1.16) determines the effect of incoherent scattering on the change of pd (in the general 
case, single and multiple scattering). 

As stated above, equation (1.16) is not applicable to the description of the process of proton 
(deuteron) transmission through the target containing light nuclei (protons (deuterons)). To describe 
multiple scattering in this case, a more general equation (1.14) should be solved. 

Further we shall concern ourselves with the deuteron passage through a carbon target, and so we 
shall use equation (1.16). 

For the sake of concreteness, let us consider the process of deuteron passage through the target 
with spinless nuclei. In this case the density matrix p(k), as well as the amplitude F(k,k') contain 
only spin variables of the scattered beam: p{k) = Pd{k). For the amplitude F{k, k'), we shall introduce 

the notation F{k, k') = f{k, k'), where f{k, k') = f{k, k'; So). 
As a result, equation (1.16) can be written as follows: 



dpd ^ T[ijy 
dz k 



if{k, k) + r{k, k)),pd{k)] + I {f{k, k) - r{k, k)),pd{k)} + 

+N j dnLfik,k')pd{k')rik',k), 

where z = vt (v is the particle velocity) is the distance traveled by the incident particle in matter. 
Hereinafter, the subscript d of the density matrix will be dropped. 

2 Scattering Amplitude 

For particles with spin 1 (deuterons), the amplitude f(k, k') can be expressed in terms of the deuteron 
spin operator S, quadrupolarization tensor Q^k and the combination of vectors k and k': 



f{k, k') = Ai + B{Sv) + CiQikfiiHk + C24fc/^H/^ife, (2-1) 

where A, B, Ci and C2 are the parameters depending on ^, z7 = [fc x k']/\ [fc x | , /T = (fc — — P| , 

fli = {k+k')/\k+k'\, the components of the tensor Q^fe are defined as: Qifc = § {SiS^ + S^Si — ^5ikl^ , 

/ is the 3x3 identity matrix. 

Write the exphcit form of the zero-angle scattering ampUtude f{k,k). In view of (2.1), we have: 

/(fc,fe) = /o(0) + /i(0)(5n)2, (2.2) 

where n = k/k is the unit vector in the direction k and the following notations are introduced: 
/o = ^ - 2Ci - 2C2, /i = 3C2. 

In the general case /o and fi arc the complex functions, and according to the optical theorem, 
the imaginary parts of /o and /i can be expressed in terms of the corresponding total cross-sections: 

Im/o(0) = A al„ lm/i(0) = |- [at] - a^,] , (2.3) 

where cr^^^, cr^J are the total scattering cross-sections for the initial spin state of the deuteron with 
a magnetic quantum number M = and M = ±1, respectively (the quantization axis z is directed 
along ft). 

The dcTitcron interacts with the target nuclei via nuclear and Coulomb interactions. The amplitude 
f{k, k') of the deuteron scattering by a target nuclei in this case can be represented in the form: 

/(A;, A; ) = fccml{k^ k) + fnucl,coul{k, k ), (2-4) 

where fcoul is the amplitude of the Coulomb scattering of the deuteron by a nucleus in the absence of 
nuclear interaction, fnucl,coul is the amplitude of scattering of Coulomb-distorted waves by a nuclear 
potential. 

The matrix elements of the amplitude /^a in the case of scattering by a fixed center are related to 
the elements of the operator T as [17]: 

777 

ha = -^y%a. (2.5) 

In the case of two interactions, the matrix elements of the operator T are defined in a standard 
manner 

Tba = {HVcoul + Vnucl\^t). (2.7) 

$a(5) describes the initial (final) state of the system "particle-nucleus" in the area, where the 
interaction is absent, the wave function if)'^ satisfies the integral equation il)'^ = $„-!- (£„ — + 
ir})~^{yccnii + ynucl)'4^t ^ where K^, is the kinetic energy operator of the incident particle (deuteron). 

Upon introducing a wave function 99^, which describes a converging wave in scattering by a 
Coulomb potential alone and corresponds to a final state ^h'- — + {^a — Kd — i'r])~^Vcoul^b ^ 
matrix elements of operator T can be represented as a sum of two terms [18]: 

ba = 'ba + >ba ' i^-^i 

where the matrix elements 7^^°"' correspond to the Coulomb interaction 7^^°"' = {^b\ Vccml Ifi): the 
part T^"'^ = <((^~| Vnud\ipt) determines the amplitude of scattering by a nuclear potential of 
waves that have been scattered by the potential Vcoui- 

Operator %iucl,coul can be expressed in terms of the operators of the Coulomb scattering and 
nuclear scattering as the following infinite series: 

Tnucl,coul Tiiucl ~l~ TnuclGoTcoul ~l~ TcoulCr o'Tnucl'^ 

+TcoulGoTnudGoTcoul + TnuclGoTcoulGoTnud + TnudGoTcoulGoTnudGoTcoul^ (2-9) 
+%oulGoTnudGoTcoulGoTnud + 



where Go is the stationary Green's function. By definition Gq = . 

Ea- Kd + 11] 

In view of the definition (2.5), equation (2.9) can be written, using the corresponding scattering 
amphtudes: 

f , jt p\ _ f ,(t p\ __ ^ f fnudjk, k")fcoui{k", fcP ,3p/ 

(27r)^m J Ek- Ek" + ^r? 

1 I fcouljk, k") fnudjk" , k') ^j'ji I 

(27r)2m J Sk- Ek" + iri 

I 1 ff fnudik- k")Ljk". P'')JnuJk'\ ^') ^jy,^jyn , 
{2iTYvn? J J {Ek - Ek" + iv)i£k - £k"' + iv) 

_l 1 f f fcoui {k,k )fnud{k ,k ) fcoul {k'" ,k') ^ k"d^ k'" — 

(27r)4m2 J J {Ek - Ek'" + iri){Ek - £fe"' + irj) 

3 Polarization characteristics of the deuterons registered by the 
detector during the beam's passage through a thin target 

3.1 Scattering in a Thin Target 

Due to a long-range character of Coulomb interaction, Coulomb scattering of deuterons by target 
nuclei occurs at small angles and the Coulomb amplitTidc for ^? <C 1 is much larger than the nuclear 
scattering amplitude. In this case, for solving equation (1.18), one can apply the perturbation theory 
[11] using as a zero approximation the solution of kinetic equation (1.18), where the collision term is 
determined only by the Coulomb interaction between the incident particle and the nuclei of matter. 
The contribution to the evolution of the spin density matrix p{k) coming from nuclear scattering 
and Coulomb-nuclear interference is considered as a correction. It should be noted that taking ac- 
count of nuclear and interference factors as perturbations is valid only for such target thicknesses z 
for which multiple nuclear scattering can be neglected [10]. In this part, in order to illustrate the 
principal patterns of relationships and simplify the form of obtained relations, we shall further an- 
alyze the characteristics of a deuteron beam for the case of a very thin target, where the change in 
the deuteron polarization state occurs only due to single scattering events in matter in addition to 
coherent scattering. 

Let us consider the process of particle transmission through the target, whose thickness z is much 
smaller than the mean free path of the deuteron in matter, i.e., z < 1/Na, a is the total cross-section 
of the deuteron scattering by a nucleus. As a consequence, in the first order perturbation theory, the 
solution of (1.19) can be represented in a form: 



p{k, z) = p{k, 0) + -N [{f{k, k) + r{k, k)), p{k, o)j z+jn[ {f{k, k) - r{k, k)),p{k, o)| z+ 

+Nf{k,ko)mrih,k)z, 

(3.1) 

where p{k, 0) is the density matrix of the beam when it enters the target, i.e, when z = 0. It describes 
the distribution over the momenta of the particles entering the target relative to the direction ko- 

In obtaining relation (3.1) it was also assumed that the initial angular distribution of the beam is 
much smaller that the characteristic angular width of the differential scattering cross-section. In this 
case the amplitude f{k, k') in the term J dQpf{k, k')p{k' ,0)f'^{k' , k) can be removed from the inte- 
grand at point k' = ko. As a result, we have J dQpf{k,k')p{k',0)f^{k',k) ~ f{k,kQ)p{0)f^{kQ,k), 
where p{0) = J dQpp{k' ,0) is the spin part of the beam's density matrix p{k,0). The term 
N J dn^,f{k,k')p{k',0)f+{k',k) or N f{k,ko)p{0)f+ {ko,k) is the contribution to the evolution of the 
density matrix, which describes single scattering of particles in the direction of k. It is known, in 
particular, that Sp/(A;, A;o)/3(0)/+(A;o, fc) is the probability for a particle to undergo a single elastic 



collision with a nucleus and get displaced by the angle corresponding to the momentum direction k 
(Sp is taking the trace over the spin variables of the deuteron) . 

Let us also consider the fact that scattering at high energies chiefly occurs at small angles ^ <C 1. 
The analysis shows that in this case, the terms in the amplitude (2.1), which are proportional to B 
and Ci, lead to insignificant depolarization of the detected beam [10] and will be dropped hereinafter. 
As a result, the amplitude / in (3.1) has a from: 

fik,k') = fo{e) + hie)iSn)\ (3.2) 

where n = k/k \s the unit vector in the direction oi k, fo = A — 2Ci — 2C2, fi = 36*2. 

Using the solution (3.1) and the explicit form of the spin structure of the amplitude f{k, k') (3.2), 
it is possible to find the dependence of the intensity and the polarization characteristics of the beam 
on the direction of the particle scattering and on the distance z traveled by the deuteron in matter. 
In a real experiment, the scattered particles are registered within a certain interval of finite momenta 
because the collimator of the detector has a finite angular width. We shall therefore consider further 
in this paper the characteristics of the beam transmitted through the target in the interval of solid 
angles Af] with respect to the initial direction of the beam propagation. In fact, due to the axial 
symmetry of the collimator, Afi is determined by the angular width of the detector collimator 2'{}det 
(for further calculations it is also assumed that "Qdet is much larger than the initial angular distribution 
of the beam). 
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Fig. 1: Scheme of scattered beam detection by the central detector. 

In the case of deuterons (particles with spin 1), the polarization state is characterized by the polar- 
ization vector P{k) = Spp{k)S and the quadrupolarization tensor Q, whose components are defined 
as Qik{k) = Spp{k)Qik- The spin density matrix p{k) can be written in the following general form: 

p{k) = h{k)i + ip(fc)l + ^Qikik)Qik, (3.3) 

where I{k) = Spp(A;). Note that alongside with the quantities P and Q, normalized spin characteristics 
of the beam are also used to describe the beam polarization. 

Denote the intensity, the polarization vector, and the quadrupolarization tensor of the beam in 
the area occupied by the detector with angular width AQ by X = J^^^ dQI{k, z), V = J^^ d^lP{k, z), 

Q = /^j^ d$7Q(A;, z), respectively. On the basis of the solution (3.1) of kinetic equation (1.18), using 
the explicit form of the density matrix (3.3) and the scattering amplitude (3.2), one can obtain the 
expression for the integral characteristics of the deuteron beam: 



l{z) =h + N {-alt + j^^ l/oPdl^) 



N 



J An J An 



2 1 

2^0 + ^{Qon)n 



V{z) = Po + iV ( -Tit + / l/ordf] )zP,+N[ --{atot - (T^tot) + 

J An J \ J An 



1 



,±1 ^0 



|A|^dJ7)n(Pon)z-^^Re 



n X (Qon)]2:, 



2;+ 



Q(z) = Qo + iV {-alt + / l/oPdf^) zQo + N {-{af^ - al^) + 2 / Re(/o/r)dQ) [ Qo+ 
V J An J \ J An J '- 



^(3n (g) n - I)/o - ^ ((Qo^) ® n + fi (g) (Qo"^)) + ^I(Qo"^)n 



z + 



l/ii'dJ] 



-(3n (g) n - I)/o- 



^ ((Qo") n + n (g) (Qon)) + +^n''Qon'' +n^n (Qon)n - ^I(Qon)n + ^Qo 
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3 27rA^ 



2 /t 



Re 



/i(0) 
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[n X Pq] (g re + n (g [n X Pq] ) z 



(3.4) 

where Pq and Qo are the polarization vector and the quadrupolarization tensor of the deuteron beam 
at entering the target, respectively: Pq = / dOP(/c,0), Qo = / dr2Q(A;,0); re = /cq/Zco is the unit 
vector in the direction of the deuteron momentum (the quantization axis) , (g is the dyadic product 
of vectors: (reig re)jj = ni-rij, is the tensor dual to vector re: {n^)ij = Eijknk, (Qo^) is the vector 
having the components (Qo^)« = Qoik^k- 

Now let us analyze the obtained solutions. First of all, pay attention to the fact that according 
to (3.4), the intensity and the polarization characteristics of the beam depend on the magnitude of 
the interval of the solid angle. As is seen, at A^l — t- 0, the contribution due to single scattering 
disappears in expressions for I, V, and Q. 

Let us give a more detailed treatment of the change in the polarization vector V depending on 
the target thickness z. Let Pq be an arbitrary, not equal to 7r/2, angle with the direction re. Choose 
a coordinate system so that the z-axis in it coincides with the direction re of the deuteron incidence 
onto the target, while the axes x and y are located in such a way that the initial polarization vector 
Pq lies in the xz plane. 




Fig. 2: Coordinate frame 



In this case Qoj/z = 0, hence the polarization vector Qore also lies in the plane defined by vectors 
Pq and re. Then from the expression for V follows that the first three terms in it describe the change 
in the component of the polarization vector, which lies in the plane (Po,re), i.e., these terms lead 
to the change in the magnitude of the polarization vector (depolarization). The last term, which is 
defined by vector [rex (Qore)], leads to the appearance in the initial polarization vector of a component 



perpendicular to the plane {Po,n), which grows with growing z, i.e., this term describes the rotation 
of vector V with respect to the direction of n. The rotation angle (peff is defined as follows: 



'Peff 



2'kN 



Re 



27riV, 



z. 



(3.5) 



As it was shown in [1, 2], the part of ^eff-, equal to ipo = — - — Re/i(0)2, is the rotation angle of 
the polarization vector due to coherent scattering of the deuterons in unpolarized matter. 

When incoherent scattering by target nuclei is taken into account, then according to (3.5), the 
rotation angle includes the additional contribution, which depends on AO. 

To provide a more detailed treatment of this dependence, it would be recalled that the deuteron 
interacts with target nuclei via nuclear and Coulomb interactions (2.4). The general spin structure of 
the corresponding amplitudes can be represented as in (3.2). This enables one to write the spinless 
and spin parts of the amplitude /, which appear in (2.4) as follows: 



U9) = a{e) + d{e), f,{9) = di{e), (3.6) 

where a{9) is the amplitude of the Coulomb scattering of the deuteron by a nucleus with the charge 
Z (the amplitude in fi (6) describing the spin part of the Coulomb interaction is small in comparison 
with di{6) [19], and so it is assumed that fcoui{9) = a{6)T); d{6) is the spin-independent and di{9) 
is the spin dependent parts of the modified nuclear amplitude. The explicit form of the dependence 
of d on ^ in the limit of small scattering angles for the case of interaction of structure particles was 
obtained, for example, in [20]. To estimate the major parameters, let us consider the most simple 
form of the dependence of d and di on 9, namely: 

d{9) = d(0)e-'='^^^'/^ d,{9) = di(0)e-*^'«2^V4^ (3.7) 
where Rd is the deuteron radius. 



3.2 Contributions to the rotation angle due to incoherent scattering 

Let us consider the terms in (3.4) which are integrated over the angular width Afi of the detector. 
According to (3.6), the functions fo{9) and fi{9) included in the obtained solutions contain the contri- 
butions from the Coulomb and nuclear amplitudes, which, in turn, have different angular dependence. 
As a result, in analyzing (3.4) for the case of high energies, we shall select such limits of -i^det that are 
bounded by the diffraction angles 9c and 9n of scattering of fast deuterons, for Coulomb and nuclear 
scattering, respectively. For a screened Coulomb potential, 9c ~ 1/kRc, where Rc is the shielding 
radius (in the case of scattering by a carbon target, i?c = 3 • 10~^ cm). When deuterons are scattered 
by nuclei whose radius is smaller than that of the deuterons, 9n ^ l/kRd- 

In view of the representation of the amplitude f{9) in (3.6), we obtain the following general 
expression for the rotation angle of the polarization vector (3.5): 



di{^)-l^ a*{9)di{9)dn- ^ d* {9)di{9)dn 
27r 27r 



z. (3.8) 



where the dependence of the nuclear amplitude on the scattering angle for high energies is determined 
by (3.7). For further analysis of the magnitude of the effect under consideration, it is convenient to 
represent ^Peff in the form 

feffi'&det) = <^0 + Vnci'&det) + fnni^det), (3-9) 

2ttN 

where ipo = — ; — Redi(0)z is the contribution of the coherent scattering to the rotation angle of the 
k 

polarization vector; (pnci'&det) = N I lm.[a*{9)di{9)\d.^lz is the part of the rotation angle, which 

J An 



describes the effect of the Coulomb-nuclear interference; ^Pnni'&det) = N I Im [d* {9)di{9)\ dilz is the 

v/An 

correction to the deuteron spin rotation from nuclear scattering. 

Numerical values of these contributions obtained, for example, for the energy of 500 MeV {k = 
0.74 • 10^^ cm~^). Then in the case of scattering by a carbon target, we have Oc ~ 10~^rad, 9n ~ 

lO^^ad. 

For the imaginary and real parts of the nuclear amplitude of deuteron scattering by a carbon 
nucleus at zero angle, take the values calculated in the eikonal approximation in the energy region 
E = 0.5 GeV: /]S'fo(0) = -5-97 + 74.81i fm, /]S'=ii(0) = -5.84 + 73.9i fm [21]. From this 
Imd = 0.75 • IQ-^^cm, Red = -0.6 • IQ-^^cu^^ j^^^ = _q qi . j^g^^ ^ q . iq-^s^.^ ^ 

result, we have ~ 1.2 • 10~*^2:. 

Let us consider the ratio ipnc/fo- Perform integration over the angular width of the detector, 
using the expression for the Coulomb amplitude in the first Born approximation. This is possible 
because for the deuteron energies considered here, in scattering by a screened Coulomb potential, the 
ratio of the real part of the amplitude in the second-order perturbation theory approximation a^^^ to 



that in the Born approximation is 
angles. 



Rea(2) /Rea(^ 



ffk kR, 



<C 1 for the whole range of scattering 



- 0.057; this estimate 



For the range of angles 9 <C 9c, the ratio is Imo/^VRea'^-* 
also remains unchanged with increasing k9Rc [22]. For k9Rc 3> 1 the ratio Ima'^^^/Rea^^-' ~ 

—2—2^ — ln{k9Rc). That is why for, e.g., 9 ^ 0^1^ tliis r3jtio is of tliG order of unity. Taking into 

n K 

account the magnitude of jlmdi/Redil ~ 7, as well as a fast decrease in the integrand due to the 
nuclear amplitude, one can demonstrate that for 9 > 9n, the contribution to ifnc/^o of the term con- 
taining lma{9) will be an order of magnitude smaller than the magnitude of (fnc/vo, calculated using 
the Coulomb amplitude in Born approximation. It should be emphasized here that this term should 
not be ignored in those rare cases, when |Imdi/Re(ii| ~ 1, as well as in the cases of calculating the 
cross-section of Coulomb-nuclear interaction and spin dichroism (for calculations in the case of high 

energies ^ 1 , it is sufficient to use the imaginary part of the Coulomb amplitude calculated in 

the second order perturbation theory). 

Thus, in the whole range of scattering angles one can obtain for tpnc/fo- 

mZe- f E. f-4 - """i"''-^ ) - Ei f-giU . (3.10) 



f°° dx 

where Ei(—z) = — — e ^ is the integral exponential function. The relative contribution to 
Jz X 

fnn/'^Q of nuclear scattering of the deuterons by the target nuclei is obtained, using the approximate 
expression (3.7): 



- \ \rJ^ - Im4 fl - e-«i'=^'^ie/2) . (3.11) 
(/^o kR\ \ Redi J V ) ^ ' 

The analysis shows that for the range of values d^et l/kRc the relationships (3.10) and (3.11) 
are much smaller than unity, i.e., the rotation angle of the polarization vector in this case is determined 
by coherent scattering of particles in matter (is determined by refraction of particles in matter). 

Note also that for stated i^det, the beam depolarization and spin dichroism will only be determined 
by the total cross-section a^^^ and a^] . 

As a result, for i^^et ^ 10~^, the system of solutions (3.4) is reduced to the solutions describing 
the evolution of the spin state of the deuteron beam only due to particle refraction in the target. 
These solutions were obtained in [1, 2]. 

With increasing 'd^^t (even for 'd^et ~ ^/kRc), singly scattered particles start affecting the po- 
larization characteristics of the beam. For instance, for i?(^et ~ 10~^ rad, the contribution of the 
Coulomb-nuclear interference to the rotation angle (pnc of the polarization vector of the deuteron 



beam is of the order of tpo, moreover, with growing i^^et ('^det ^ ^n), according to (3.10), (pnc grows 
logarithmically: (fnc/vo ~ ln(A;i?c'!?(iet)- 

3.3 Rotation angle including the Coulomb contribution to the spin part of the 
nuclear amplitude di{0) 

It should be noted that if in (3.8), integration over dfi is made over the entire solid angle (a An 
experimental geometry), then according to the analysis given in [11, 12], the considered interference 
term and the contribution to the amplitude di{0) due to distortion of the incident waves by Coulomb 
interaction compensate one another, i.e., the sum of the contributions ipo + ^nc{'&det = t^) is in this 
case determined by the spin-dependent part of a pure nuclear amplitude of scattering at zero angle 

(^Q = yjQ + ¥'nc(7r) = — —Kedi. This will be demonstrated below. 

The modified nuclear amplitude fnud,coui{k,k') is associated with "pure" Coulomb and "pure" 
nuclear amplitudes by formula (2.10). Taking into account that the solutions of (3.4) and (peff include 
only the spinless part of the scattering amplitude a{6), according to (2.9), the spin-dependent part of 
di can be written in the form: 

d,{k,k)-d,{k,k) ^2^^,^y e,-e,. + ir, ^ {2^fm J e,-e^„^i^ ^^'^^^ 

The amplitude d'^ is a "pure" nuclear spin-dependent amplitude of scattering at angle d. 

Substituting the explicit form of the amplitude di{k,k') from (3.12) into the expression for (feff 
in (3.8), and again taking into account only the first degree of the product of d'l and a, one obtains 



27r7V 



d'lio) + ^ / a{e)d[{e)dn f a*{e)d[{e)dn f d*{e)d[{e)dn 

27r J 27r J^^ 2tt J 



(3.13) 

Integration in (3.12) was carried out subject to the following representation of the stationary Green's 
function: 

^— - = P ^ iT^5{ek - £fe')> (3.14) 

where P is the principal value integration. 

As is seen from (3.13), when the detector registers particles in the 47r solid angle geometry (Afi = 
47r), the Coulomb-nuclear contributions are mutually canceled. The compensation, however, only 
occurs when we can confine ourselves to the consideration of the Coulomb scattering amplitude in 
the first Born approximation. 

Allowing for the Coulomb contribution to cii(O), equation (3.9) can be written as: 

¥'e// = v'q + Vnl + Vnn, (3.15) 

where ip^°^ is the total contribution of the Coulomb-nuclear interference to the rotation angle of the 
polarization vector. Its magnitude is determined by the sum of the second and third terms in (3.13). 
Within the limits of small scattering angles, the explicit form of ip]"^ reads as follows: 



ip^°l = 2itNz / Im [a{e)d'^{e)\ 9d9 = -2TrNlmd[—^Ei 



4m 



z. 



In this case the relative contribution of the Coulomb-nuclear interference to the rotation angle is 
determined by formula 

^iot_ lmd[mZe^^.r k'R'AA ^^^^^ 



For high energies, when nuclear scattering has a pronounced diffraction character, one can estimate 
the magnitude of the lower limit of the angle ^^XT''' which the interference contributions can be 



considered to compensate one another with the selected accuracy. We have that for the considered 

energy and the form of nuclear amplitude, the total contribution of the Coulomb-nuclear interference 
can be neglected even at i^det On (the diffraction angle of nuclear scattering) (l*/?^"*! order of 

magnitude smaller than ip^). 

For example, if-^det = 0.8- lO^^ rad, ipnc/v>o = -0-2. Prom this it can be assumed that -e^j;^" ~ 

If the collimator of the detector registers particles moving within a certain solid angle AQ corre- 
sponding to i3ciet ^ ^'di^'^j the stated Coulomb- nuclear terms will not be compensated. For example, 
when ■ffdet — 0.4 • 10"^ rad, the part of the rotation angle (^^°* is comparable in magnitude with <^q. 

It should be emphasized that the minimum value of "ddet, when the Coulomb-nuclear contributions 
are compensated (i?^°7'')) depends on the deuteron energy ('^ l/V^) as well as on the type of the 
dependence of the spin part of the nuclear amplitude on the scattering angles. 

Consider now the nuclear contribution to the rotation of the deuteron spin. The estimate of 
relation (3.11) for -^det ~ Gn gives ipnn/'Po - -10"^. With growing i^det, namely for -ddet > '^/kRd, 
the ratio \(pnn/fo\ achieves its maximum value: \<fnn/ fol — 0-3, while the magnitude of the ratio 
\^nn/^'o\ is of the order of 10~^. 

Let us mention the following possibility of experimental measuring of the spin-dependent part of 
the zero-angle nuclear amplitude. Measuring the rotation angle of the deuteron polarization vector for 
two arbitrary values of the angles i?! and 'i?2 of the detector and considering the difference (peffi'&i) — 
(peff{'&2), one obtains that ipeff{'&i) — (Peff{^2) ~ ^nc{^i) — ^nc{'&2)- Using equation (3.10), one can 
thus find the value of Imdi and, hence the value of (/?„c for any "ddet- 

Knowing them and taking into account that |(/J„„| <C \^nc\ for any ^det-, it is possible, using (3.9), to 
estimate the magnitude of the rotation angle (po due coherent scattering and, hence, the magnitude of 
Redi. This part of di(0) can also be obtained directly measuring the rotation angle of the polarization 
vector for -ddet ^ 1/kRc- Basing on the assumption that the Coulomb-nuclear contributions in ip^ff 
are compensated and choosing the angle of the detector so that ■ddet ^ "^"det^ ^ can obtain the 
magnitTidc of Red'^ . 

The domain of applicability of the solutions of (3.4) corresponds to such target thicknesses z for 
which multiple scattering in matter can be neglected, i.e., z is much smaller than the mean free path 
of the deuterons due to Coulomb scattering: z < l/Nacoul- For energies 0.5^1 GeV, thickness z is of 
the order of 10""^ cm. 

In the next section we shall demonstrate the validity of the assertion that in the cases when 
z is much larger than the mean free path l/Nacoul (i-e., in the targets where multiple Coulomb 
scattering takes place) , the contribution of the Coulomb-nuclear interference to the rotation angle of 
the polarization vector also depends on the angle of the detector collimator. 



4 Polarization characteristics of deuterons under the conditions of 
multiple Coulomb scattering in matter 

With increasing target thickness, at least the condition of the smallness of the mean free path due 
to Coulomb interaction Na^oul^ <C 1 is violated. In this case the solution of the kinetic equation, 
which describes the Coulomb scattering of a deuteron by the nuclei of matter should be found. It 
will describe the beam distribution due to single and multiple Coulomb collisions of particles in the 
target. 

As has been stated above, if in this case z ^ l/Nanud, which is realized in the majority of practical 
cases, one can find the nuclear and Coulomb nuclear contributions in the scope of the perturbation 
theory. 



4.1 Solution of the kinetic equation 



Write equation (1.18) as follows: 

dp{k) n N 



= -Nal,m - f - al,) {{Snf,p{k)} + ^iVRe/i(0) [{Snf,p{k) 



dz 



+ 



(4.1) 



where a°^^ = Sp f AVLj:,f{k,k')pQ^^f^{k,k') + cr"'='=\ /5q*^ is the deuteron density matrix describing 
the initial state of the beam for M = and M = ±1, correspondingly; ft is the unit vector in the 
direction of fc; (Jy. denotes the part of the total cross-section, which is responsible for all possible 
inelastic interactions between the deuteron and the nuclei of matter, including nuclear reactions. 

Substituting (2.4) into the expressions for the total cross-sections, the elastic part of cr°^f ^ can 
be represented as a sum of the terms describing the cross-section of the Coulomb scattering, the 
pure nuclear cross-section and the cross-section of the interference between the Coulomb and nuclear 



interactions: a 



0,±1 



tot coul ' nucl,coul ' nucl 

Use the following notations: a'^j^ = a! ^ 



(T°of ^ in the form: 



, 0,±1 _|_ 0,: 

nucl, coul '^nucl ~^ 



o,±i 
'tot 



^0,±1 , J,±l 



'NC ' 



Thus, represent the cross-sections 



(4.2) 



where a°„t) = Sp J dn^,J^^ul{lk')pr' fLli^' a%%' = Sp J d^^,, {fcoul{kApTfLcl,c<^li^'.k) 



+ 



+fnucl,coul(k,k')pQ^ f^gy^lik' ,k) + fnucl,coul(k,k')pQ'^ fnucLcoul^^' ^ ) 



o,±i 
r 



In view of the notations introduced above and the representation of the scattering amplitude (2.4), 
rewrite kinetic equation (4.1) as: 



-Nal^^iPik) - ^ {a^i - {iSnf,p{k)} + ^NRefr\0) [{Snf,pik) 



dp{k) 
dz 

+N I dn^,f^i{k, k')p{k')f+Jk', k)- 
N 



+ 



+ 



(4.3) 



-Na^cm - y Kc - ^%c) {iSn?,p{k)} + ^iVRe/r^(0) [{Sn)\p{k) 
+N J d% {houi{lk')p{ky^:^^^^^^^^^^^ 

+N J dJ^g,/w,co«K^,^')p(^')/:„d,co«K^''^)' 

As stated above, the solution of equation (4.3) can be presented as follows: 

p(^,z) = p(°)(^,z)+pW(^,z) + ..., (4.4) 

where p^^\k, z) is the zero approximation, which is the solution of kinetic equation (4.3) and describes 
only the Coulomb interaction between the deuteron beam and the nuclei of matter; p^^^ {k, z) is the 
first order perturbation theory correction allowing for the contribution of nuclear scattering to the 
evolution of the polarization characteristics of the beam. 
The equation for p^^^ (z) has a form: 



+N J dn^Jcoul{k,k')p^''\k')f:^i{k,k'), 



27ri 



{iSnf,p^'\k)} + ^^ArRe/r'(0) {Snf^p^'Hk) 



(4.5) 



where the spin structure fcoul{k,k') in the general case has a form (2.1). 

It can be shown [10] that the terms in the Coulomb amplitude (2.1), which are proportional 
to B, Ci, and to the part of C2 depending on 9, lead to depolarization of the registered beam. 

The magnitude of such depolarization is determined by bg9'^z {hg = ^ ^ ^ - ^ + — g is the 

/dcT ^ ^ 

9^—^^dQ is the average squared angle of Coulomb scattering). With 

due account of the small values of the scattering angles due to Coulomb scattering at high energies, 
this quantity will be insignificant as compared with the contribution from the scalar part of the 
amplitude (2.1). 

For this reason, in the case of Coulomb interaction we can confine ourselves to considering the 
scattering amplitude of the form (3.2). 

Moreover, the term in (3.2), which is equal to /f°"'(0)(S'n)2 and describes a coherent rotation of 
the polarization vector of the deuteron beam due to the elastic Coulomb scattering at zero angle, is 
also small and thus can be neglected [19]. 

As a result, equation (4.5) is rewritten as follows: 



dz 



Nal,lP^''\k,z) + N I dn^, \a{k,k')\%^''\k',z), (4.6) 



where a{k,k') denotes the spinless part of the amplitude fcoul{k,k'). 

Expression (4.6) is an integro-differential equation. Its solution in the limit of small scattering an- 
gles can be obtained by expanding the function (k') over a relatively small parameter q (transferred 
momentum) [23]: 



p(^){k',z)^p^"){k,z) + ^q,+ 



+ ^ 



2 dki 



+ ^ 



1 ^2^(0) 



2 aA;2 



(4.7) 



where Qx = qcostp, qy = qsincp, 2± k9, where 9 is the scattering angle (the angle between vector k 
and the z-axis). 

Substituting the expansion (4.7) into (4.6) and integrating over the azimuth angle cp, one can 
obtain 



dp(^\k,z) 02 r ^ ^^^A-^mwr X 

= -^^P^^'HKz) + — AApW(A;,z) + 

^2 92 

operator A acts on the transverse components of vector k: A = 7^—^ + 77-^, where ft 



,2d<ycoul_^^^ 04 = / 9^^^dn, and etc. 



9n2 9n2 ' 



(4.8) 

k/k; 



9"^ = N 

' " dn " " y " dn 

If only the first term on the right-hand side of equation (4.8) is taken into account, the integro- 
differential equation (4.6) is reduced to a parabolic differential equation. The limiting angle for such 
approximation is obtained from the condition [23]: 



9 

''max ^ 



The solution of equation 



dpW(fc,, 

dz 



16ApW(fc,, 
AAp(0)(fc,; 



^ApW(^,.) 



for the initial condition p^^\k,z = 0) = PQS{nx)S{ny) and an infinite medium is 

p^^\k,z) = posgc{k,z), 



(4.9) 



(4.10) 



(4.11) 



where 111 

Pos = -^IqI + 2^oS + -QoikQik 

(see equation (3.3)), gc{k,z) = e e'^z . The unit vector n is counted from an arbitrary 

direction ng (then direct the z-axis along no), {n — no)'^ = 9'^, where 9 is the angle between k and 
the z-axis directed along no; 9"^ is the average squared angle of scattering per unit path length. The 
average squared angle of scattering at depth z {9'^z) will further be denoted by 9'^, and its square 
root, by 9z- 

The applicability condition (4.9) of this solution is equivalent to the inequality 9 <^9z- It should 
also be stated that for large 9, the solution of equation (4.6) will mainly behave as ~ [24]. Thus, 
in the approximate equation (4.10), scattering at large angles in a single scattering event is ignored. 
Indeed, for 9 » 9^, the solution (4.11) decreases exponentially, while the next term, corresponding 
to the solution of the initial equation (4.6), decreases according to a power law 9~^. The angular 
distribution (4.11), therefore, does not describe the characteristics of particles scattered at angles 
9:$>9z. 

The expression for the density matrix in the zero-order approximation (4.11) enables obtaining 
the correction in the first order perturbation theory p^^\k,z): 



{k, ^)=n£ J G{k -k";z- z') {-a^cp"-'^ {k" , z') - \ {a% - a%c) {Uk", k"), {k" , z') } 
+ ^iVRe/r^'(0) [/,(fe",fc"),p(°)(^",^')]) dn^„dz'+ 
+N j G{k-k"-z-z') {hUk'',k')p('\k',z')f:^^i^^^^^^^ 



+ 



+ fnuci,coui{k", k')p^^\k', z')f+{k', k")) dnpdn,„dz' 



+N y" J G{k - k"; z - z')fnucl,coul{k", k')p^''^ {k', z')f+^^^^i{k', k")dnpdnj^„dz', 

(4.12) 

^ ^ 1 (n-n")^ 

where G{k—k"; z—z') is the Green function of equation (4.10): G{k—k"; z—z') = e e2|z-z'| 

Tr9'^\z — z'\ 

The amplitudes of Coulomb and nuclear scattering in the general case are expressed by formula (2.1). 

Thus, the solution of kinetic equation (4.1) allowing for nuclear interaction has the following 
general form: 

p{k,z)c^p^''\k,z)+p^^\k,z) , (4.13) 

where the main contribution p^^^k, z) is defined by expression (4.11), the correction in the first order 
perturbation theory p^^\k,z) equals (4.12). 

To represent the solution in the explicit form, we shall confine ourselves to considering small 
scattering angles. In this case one can obtain the following approximate expression for the amplitude 
(2.1): 

f{k, k') = A{x) 1 + bx{Sv) + cix^Qikf^iPk + C2(l - X^/4^)Qikf^iil^ik , (4.14) 

where the angular dependence of functions b{x), ci(x), and C2(x) was supposed to be b{x) ~ sinx, 
ci (x) ~ sin^ X, C2 (x) ~ cos^ x/2. 

Integration of the right- and left-hand sides of the solution of (4.13) over the final momentum k 
in the limits corresponding to the variation of the axial angle from to 2tt and the polar angle, from 
to a certain 'd^et gives the expression for the density matrix of the beam that has passed through 
the area occupied by the detector with angular width AO. 

As has been stated above, for high energy deuterons, elastic scattering occurs mainly at small 
angles. In this case in the solution of (4.13) the general structure of the scattering amplitude in 
the form (4.14) should be used for the nuclear amplitude. It can be shown, however, that under 



the assumption of axial symmetry of the conimator of the detector, the correction terms in the 
polarization characteristics can be neglected with high accuracy. These polarization characteristics 
correspond to the terms in (4.14), which are proportional to x &iid ■ This assertion is also valid 
for such target thicknesses and deuteron energies, for which the following inequalities are fulfilled: 
< 1, l/kRd < 1, l/kRc < 1. 

As a consequence, in the expression for the density matrix (4.13), the amplitude fnucl,coul{k,k') 
of scattering due to nuclear interaction is taken equal to: 



fnud,coul{k, k') = d{e) + di{e){Sno) , 



(4.15) 



no is the unit vector directed along the 2;-axis, d{9) is the spin-independent and d\{9) is the spin- 
dependent parts of the nuclear amplitude, 9 is the angle between k and k' . 

For the Coulomb scattering amplitude fcoul{k,k'), its spinless part a{9) can be used with high 
accuracy. 

Integral characteristics of the deuteron beam are obtained by substituting the explicit form of 
the density matrix (3.3) and the scattering amplitude (4.15) into the solution (4.13). As a result, we 
have: 



I{z) = (1 - e «i )Io + ^ilo + (6 + 6 



2 1^ 



"det 
fl2 



_ I r _ _-| _ 2 

P{z) = (1 - e «i )Po + ^iPo + [i=b + no(Pono)J + 6^o(J'ono)z - -<^e//[no x (Qono)], 



get 
fl2 



Q{z) = {l-e «i)Qo + eiQo+ 
+6 



Qo + ^(3no <8) no - I)i"o - ^ ((Qo^o) <8) no + no <X) (Qo^o)) + ^I(Qo«o)no 



+ 



+6 



^(3no (g) no - I)/o- ^ ((Qono) no + no (Qono)) + ^noQono + 



+ no (8) no (Qono)no - ^I(Qono)no + ^Qo 



-^ipeff [[no X Po] <8) no + no (8) [no x Po]j ■ 



(4.16) 



The parameters ^i, ^2, and ^3 are defined in terms of the relations: 

z 

An 

G{k - k"; z - z') (a{k", k')d+{k', k") + d{k" , k')a+{k', P')) 9c{k' , z )dil^,M^Mj:dz'+ 
+N J J J J G{k-k'';z- z')d{k'',k')d+{k',k'')gc{k',z')dnp,d^pdn^dz', 



An 

AH 



AH 



z 

U^) = -N{a%'c - <^Nc) J J / ^(^ - ^ - ^')9c{k", z')dn'idnj^dz'+ 

An 

z 

+^ J J J J Gik-k";z- z') (^di{k\k')a+{k\k'') + a{k\k')d+{k\k'')^ Qc^^^ 

AH 

z 

+^ J J J J ^ ~ (^(^"' ^')dtik', k") + di{k", k')d+{k', k")^ g^k', z')dn'idn^,d^gdz' , 



An 



Uz) = N j j j J G{k-k";z- z')diik", k')dtik', k")gc{k', z')dVL^,4^^4^^dz' . 
AH 

(4.17) 

The general structure of the solution (4.16) is written in the same way as the solution (3.4) for a 
thin target. For this reason, all the conclusions derived in analyzing each of the terms appearing in 
the expression for the polarization vector in (3.4) are automatically (naturally) extended to the case 
of a thick target. 

The rotation angle of the polarization vector with respect to vector no appears equal to 

2 

'^e// = ^ y" j j G{k- k"; z - z') (di{k, k) + d+(fe, k))gc{k", z')dn^„d^^dz'- 
An 

z 

/// j ^^^-^"'^^-^') [di{k'',k')a+{k',k'')-a{k'',k')d+{k',k'')^gc{k^ 



AQ 

z 



/ / / \ ^ ~ ^'^ {di{k"^ k')d+{k', k") + d{k", k')d+{k', fc")) 9c{k', z')dn^„d^j:,d^^dz' . 

AQ 

(4.18) 

4.2 Analysis of the problem of compensation of the Coulomb-nuclear contribu- 
tions in the expression for the rotation angle of the poleirization vector 

Let us consider in detail the contributions that are included in the expression for the rotation angle 
of the polarization vector of deuterons. The expression between the first brackets that contain the 
nuclear amplitude di is the part of the rotation angle, which is due to the coherent scattering of 
deuterons in matter; the expression between the second brackets, which consists of the product of di 
and the Coulomb amplitude a, describes the effect of the Coulomb-nuclear interference; the last part 
in (^e//) which contains the prodiict of the spin-dependent and spinless parts of the nuclear amplitude, 
is the correction to the deuteron spin rotation from nuclear scattering. 



It would be recalled, however, that nuclear amplitudes d and di are actually the scattering ampli- 
tudes modified by the Coulomb interaction. The relation of fnucl,coul{k,k') to "pure" Coulomb and 
"pure" nuclear amplitudes is expressed by a general (2.9) or an approximate (3.12) formula. 

Substitution of di{k,k') from (3.12) into the expression for <^e// (4-18) in the accepted approxi- 
mation gives 



^eff = ^Re4(0) / / I G{k- k"; z - z')gc{k", z')dn^„dn^dz'+ 



k 

Af2 



iN 

An 




G{k - k"; z - z') (d[{k, k')a{k', k) + a{k, k')d[{k', k) - 
-d[*{k, k')a*{k', k) - a*{k, k')d[*{k', fc)) g^k", z')dn'^jnpdn^dz' - 

Gik - %"■ z - z') (d\{k", k')a+{k', k") - a{k", k')d+{k', fc")) gdk', z')dQ.^,M'^,d^j^dz' - 




_iN 
2~ 

An 

z 

/ / / j ^ ~ l^')d'^(.k', k") + d'{k", k')d'i{k', fc")) gc{k', z')dQ^„d^^4^^dz' . 

An 

(4.19) 

Expand the Green function G{k — k'; z — z') and the function gc{k,z') into a series of orthogonal 
Legendre polynomials: 

G{k -k';z- z') = 5Z + ^) ^"(-^ - '2')^'n(cos^?), 

5c(P, ^) = ^ + -j9cK^)^'Kcos t?'), 

where •& is the angle between vectors k and k'\ d' is the angle between vector k' and the z-axis. 

In substitution of relations (4.20) into (4.19), one should take into account the condition of the 
orthogonality of the polynomials: 

I P-(cos^)Pr(cos^) sin^d^ = -^Ij^-^.r. (4.21) 



and summation theorem: 

P„(cosx) = Pn(cos^)P„(cost?i) + 2 V p-^P™(cos^?)P™(cosi?i)c (4.22) 

m=l ^ ' 

where % is the angle between the two vectors that make with the z-axis the angles d and ■i?!, respec- 
tively. Considering small deviation angles of particles from the initial direction (the z-axis), let us 
take account of the fact that large n and I play the main part in the expansions (4.20). So in view of 
of the below relation, one can pass from the expansion into a series of Legendre functions to that of 
Bessel functions 

lim P„(cos^) = Jo(i?), (4.23) 

where Jo(t?) is the zero-order Bessel function. In this approximation the expansion coefficients Gn{z — 
z') and gcl{z) have a form: 

GJz-z') = -e-'^^^^-^'\ 



As a result, expression (4.19) can be rewritten as follows: 



AH 

/ / ('^i^^' k) + a{k, k')d[{k', k) - d'{{k, k')a*{k', k) - a*{k, k')d'^{k', k)^gc{k, z)dnpdn^- 

AQ 

--Y^ j j ('^i^^' ^>*^^^ - «(^' ^')d'^{k, k')^gc{k', z)d^^,d^^- 

AQ 
AH 

(4.25) 

Let now integration over the directions of the final momentum k be made in the domain of the 
entire range of variation in the scattering angle. In this case in (4.25), it is possible to replace the 
integration variables k by k' and vice versa. 

For such deuteron energies at which only Coulomb amplitude in the first Born approximation can 
be taken into account, the Coulomb- nuclear contributions in the expression for the rotation angle 
compensate one another. 

As a result, registration of the scattered and transmitted particles in a 47r experimental geometry 
gives the following magnitude of ipef / : 



2'kN ^ 



d'M d[{k,k')d'*{k,k')gcik',z)dn^,dn^ 



(4.26) 



In the case under consideration, the rotation angle is determined by the sum of two terms: the 
coherent contribution, depending only on a pure nuclear amplitude of scattering at zero angle and 
the incoherent contribution due to single nuclear scattering. 

In a real experiment, however, this case of a 47r-gcometry is generally not realized. That is why, 
integration over k in (4.17) and (4.18) or (4.19) should be made within finite limits. 

4.3 Calculation of the peirameters of integral chciracteristics of the beam. 

Integration over the angular variables in the parameters ^i, ^3, and V'e// is performed using the 
expansion of the Green function G{k — k';z — z') and the function gc{k, z) into a series of orthogonal 
Legendre polynomials according to formulas (4.20)-(4.24). 

In the explicit form, the parameters of the system (4.16) are as follows: 

^2 OO 

Uz) = -Nalfc{^-e'^)z + 27rNz J P(x; ^?det,"^) ( 2Re[a(x)(i*(x)] + \d{x)\^)xdx, 



,o2 CX) 

Uz) = -Nia%'c - a%c){^ - e ^^)z + 2'kNz j P{x; ^det, ( Re[a(x)dt (x)] + 



+ Re[d{x)dUx)])xdx, 

OO 

Uz) = 27rNz J P{x;^det,'0l)\di{x)Wdx, 



(4.27) 



27riV 



l_e"#j -ik J P{x;^det,'0!){a{x)di{x) + d*{x)di{x)) xdx 







(4.28) 



where the integral P{x;'&det, G^) = 



2 

'0d'&ndne~^^^Jo{n'd)Jo{nx) is denoted in terms of the 







function P{x;'&det, ^z) The last expression is integrated over n. As a result, the function P can 
be represented as the integral over ■& of the expression containing a modified zero-order Bessel 



function {Io{z) = Joiiz)): Piw'^det 



"i Iq [ 1 ^?o?7? . The explicit form 



of the introduced function, which can be written in a form of the following infinite series is 

m+l 



m=0 ^ ' 

sian hypergeometric function. 



m I det 



-m; —1 — m; 1; 



where F is the gaus- 



4.4 Analysis of the contributions to the rotation angle yje// for the case of scat- 
tering by a thick target 

Rewrite expression (4.28) for two limiting values of the angle ddet of the collimator of the detector: 
•^det < Gz and t?det » Gz- 

Before considering the angles -iJ^ct ^ let us indicate the following points. As has been stated 
above, the solution (4.11), which describes the angular distribution of the deuterons due to Coulomb 
interaction holds true for the scattering angles 9 <^ 6z- Since we analyze the integral characteristics 
of the beam, the solution (4.11) can also be used for i^det S> Oz- The stated approximation means 
that the contribution of singly scattered particles due to Coulomb interaction is neglected. Moreover, 

1.. ..... . e .a.e »..*«Uon . .... / .-(M).. - a .... eo..W«„„ . 



^det 

p^^\k,z)(Kl (the major contribution to the integral comes from small 9 ). 



In the first case, in integration of the second term in (4.8), it can be assumed that Jo(n^9) k. 1. 
Let ^\ff denote the rotation angle of the beam that has passed through the area occupied by the 
detector with angular width i?det 9z- The explicit expression for ^\ff has a form: 



/ 



^SiRe 
k 91 



dM-ik J xdxe'^ {a*{x)diix)+d*{x)diix)) 



(4.29) 



For another limiting case -^det 3> ^z, relation (4.28) is written as: 

27rA^, 



k 



-Re 



di{0)-ik / Xdx{a*{x)di{x) + d*ix)d,{x) 



z. 



(4.30) 



One can easily see that in the case when the whole beam gets into the detector, the expression for 
the rotation angle ip^^j^j: is similar to that for (feff, which was obtained in describing the passage of 
the deuteron beam through a very thin target in the case of a wide experimental geometry (3.8). 
As it has been done in the previous section, let us represent ipeff as a sum of three terms: 



^eff{-&det) = Voi'&det) + fnc{-&det) + ^nni'&det), 



(4.31) 



where <^o('*?(iet) = — — Redi(O) ~ ^ ^det/^i^ z is the contribution of the coherent scattering to 



the rotation of the polarization vector; ifnci'&det) = 27rA^Im 



P{x;^det, Ol)a {x)di{x)xdx 



z is the 



part of the rotation angle, which describes the effect of the Coulomb-nuclear interference; ^nni^&det) 



27r7VIm 



z is the correction to the rotation of the deuteron spin from 



a pure nuclear scattering. Obtain the numerical values of these contributions for the energy of 
500 MeV, k = 0.74 • lO^^cm"^. Then we have kRd = 31.82 and the value of the average squared 

angle of deuteron scattering by a carbon target per unit length 6'^ = 16ttNZ ( — I ln(137Z ' ) 

\pvj 

equals 0.23 • lO^^cm^^. The maximum value of the rotation angle ifo is obtained at -ddet ^ Oz- 
ifo ~ 0.12 • W^^z. 

Consider the relation ipnc/vo- Calculate the explicit form of the Coulomb-nuclear correction, 
using the expression for the Coulomb amplitude corresponding to scattering by a screened Coulomb 
potential in the first Born approximation: 



a{e) = - 

In the case of scattering by carbon, kR^o-ai 
whole range of scattering angles: 



mZe^ 



-R. 



1 



^2 n + A;2i?2^2 • 



(4.32) 



2 • 10^. Using (3.7) and (4.32), one can obtain for the 



Hedi TT?k 



1 



det I 2 



1 



1 + /£2i?2^2 



'Xdx 



(4.33) 



Plot this dependence for three values of 01 (since the value of energy is fixed, this corresponds to 
three values of the target thickness z\. 



Vnc 
90 



7.0 



6.5 



6.0 




0.02 0.04 0.06 0.08 0.10 



^det 



Fig. 3: Relative contribution of the Coulomb- nuclear interference as a function of the angle of 
the detector for 9^^ = 0.23 • 10~^ cm~^ (solid curve), 6*22 = 0.23 • 10~^ cm~^ (dashed curve) and 
"^3 = 0.23 • 10-2 cm-i (dotted curve). 

A rather simple analytical expression for (pnc/^o can be obtained for two limiting values of 'ddet- 
Let t?det 0z- From (4.29) we have that for small values of the angles of the detector, the relation 
^nc/^o is the function independent of 'ffdet- 



fnc 



Imdi mZe^ 
Redi h^k 



Ei 




exp 




(4.34) 



When -i^det 00 the expression (pnc/^o tends to its limiting value (see Fig. 1): 



Redi h^k 



Ei 



4i?2 



(4.35) 



According to (4.34) for such energies and target thicknesses when 9^ ^ 9^ the relation {'^nc/'-Pof 
tends to ((/?nc/¥'o)") i-e-, weakly depends on "ddet- If the parameter 9z <C 9n-, then for -t^j^et — )• the 
relation {ipnc/^oY decreases with decreasing target thickness. The corresponding contribution can 

fi^k Rcdi\ 



be neglected in the case when 9^ <C 9^ exp C 



I , where C is the Euler-Mascheroni 



mZe^ Imdi / 

constant, C =0.5772. For the energy of 500 MeV, the value of z that satisfies this condition is 

32 „„„ ( ri 1 ^YiQ value of tpnc at "ddet ^ is of the order of 



z < 10-*^ cm. When 91 > 9^ exp C 



mZe^ Imdi 
or greater than the main contribution ipQ. 

With increasing '&det-> the relation ^nc/^o 
Imdi mZcP' 



Redi h^k 



C + In 



4i?2 



get thickness. From this fo^ 



grows, approaching (v^nc/v^o) = 
and this value of the contribution does not depend on the tar- 



flows that when t^det ^ &z-, the contribution of the Coulomb-nuclear 
scattering cannot be neglected for any z: cpnc is one order of magnitude larger than the main 
contribution (pg. 



4.5 Rotation angle including the Coulomb contribution to the spin part of nuclear 
amplitude di{0) 

Let us consider in more detail the expression for the rotation angle including the Coulomb contribution 
to the amplitude di{0). For the analysis of the total contribution of the Coulomb-nuclear interference, 
it is convenient to write equation (4.31) in a form: 



nrii 



(4.36) 



where (p °l is the total contribution of the Coulomb-nuclear interference into the rotation angle of 
the polarization vector. Within the limit of small scattering angles, its explicit form reads: : ip\ 

oo 



tot 
nc 



2ttN mZe^ 



Imo?i < 



J fix; 



1 



4 xdx + Ei 



3. 

ml 



,92 



1 + k'^Rlx^ 

Basing on the above considerations, one may immediately see that when the condition 9z S> 9n is 
fulfilled, the total Coulomb-miclcar contribution vanishes. 

Let us analyze the behavior of the Coulomb-nuclear interference when 0^ <C ^n- With this aim in 
view, consider the relative contribution of (/5^°*: 



> z. 



tot 
rnc 



+ Ei 



Imd'^ mZe^ 



Red[ h?k 



1 



+ k'^Rlx'^ 



Xdx+ 



3. 

4i?2 



(4.37) 



where Hed'i is the spin-dependent part of a pure nuclear amplitude of scattering. According to (3.12), 
the relation between the amplitudes Red'^^ and RecZi is determined by formula 



TTlZ 

Red'i = Redi + Imdi Ei 



3: 

ml 



(4.38) 



The calculated value is Redi/Red'^ = 0.12. Using this estimate, let us plot the dependence v'nc/'^o 
for three values of the average squared angle 9^ of multiple scattering at depth z: 
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Fig. 4: Relative contribution of the Coulomb-nuclear interference V'nc / '^d ^ function of the angle 
of the detector for 6^^ = 0.23 • 10""^ cm~^ (solid curve), 6*^2 = 0-23 • 10~^ cm~^ (dashed curve) and 
0.23 • 10~^ cm~^ (thick solid curve). 



It is seen in Fig. 4 that the maximum absolute value of \(Pnc \ is achieved at the angles of the 
detector i^det ^ (^z- 



tot 
^0 



Ei 



Imdi mZe^ 
Red; h'^k 



3. 

Am 



Ei 




exp 



1 



k'^Rl I ^ 4 



(4.39) 



For considered energy, at target thickness z < 0.01 cm, the magnitude of \ipnc\ is of the order 
of ipQ. With further decrease in z, the contribution of \(Pnc \ grows and reaches its maximum value 

TT^jT-, — El ^ at z — )■ 0. It IS easy to see that it is exactly equal to the maximum absolute 

Rea^ n^K [ 4x2^ J 

value of the relative value of the total Coulomb- nuclear contribution \'f*nc/v'o\ the case of a thin 
target (3.16). The maximum value of this quantity is also achieved for small 'ddet^ namely, within the 
limit ddet ~^ 0. Moreover, for a thin target, the dependence of the relation </^^°*/vo on the detector 
angle "Qdet is, in fact, the limit to which the Coulomb-nuclear contribution tends when deuterons are 
scattered by a thick target, if the parameter 61 ^ Q (dotted line in Fig. 4). 

Indeed, at 01 — 0, the introduced function /(x; i^dot, ^2) becomes equal to f{x]^det,0) = 

J 6{'d — x)d'Q. As a result, for the rotation angle '^eff defined by formula (4.28), we obtain the 


expression for V'e// in the case of deuteron scattering by a thin target (3.8). It is clear because 
within the limit 91 — >• 0, the contribution of multiple scattering can be neglected for given energies 
and (or) target thicknesses. This statement is true even for Oz <C Oc, i.e., for the case when the 
value of the average squared angle of multiple scattering is much smaller than the diffraction angle 
for Coulomb scattering. In this case, similarly to the case of scattering by a thin target, one can 
introduce the detector angle f?^°7'' at which one may consider that the interference contributions 
compensate each other. Thus, for 'd^et ^ the contribution of is comparable in magnitude to 
(^0, i.e., the Coulomb-nuclear terms are not compensated. When ^det > Gn, the total contribution of 
the Coulomb-nuclear interference can be neglected. 

On the other hand, when 6z ^ 6c, the dependence of (/9^°*/¥'o "^det becomes more smooth (Fig. 



4). In this case the the maximum absolute value of \ip^°l^j ip'^ diminishes with growing parameter ff^: 



tot 
^0 



Red; h?k 



In 




mi 



(4.40) 



Prom this follows that at the energy of 500 MeV and the target thickness z ^> 1 cm, up to the 
accuracy of 10%, one may consider that the Coulomb-nuclear interference due to incoherent scattering 
compensates the Coulomb-nuclear contribution to the amplitude di{Q) for all values of "ddet- The 
expression for the rotation angle '■Peff will also take quite a simple form (below, it is shown that the 
nuclear part of ipnn can as well be neglected in comparison with c^q): 



2tiN 



)z. 



Using (4.38), we obtain that the maximum value of this quantity is (fQ ~ 10 



-2 



z, I.e. 



(4.41) 



one order of 



magnitude larger than ipQ. 

Thus, for the energies in the region 0.1 1 GeV, one can indicate two main domains of variability 
of the value of the multiple scattering parameter 9^ ■ 

1. 6z <^ Oc- Here the dependence of the function ^p^nc/ y^'o the angle -i^^ei of the detector is close 
to the similar dependence in the case of a thin target. 

2- ^2 ^ ^c- In this range of energies and target thicknesses the dependence of the magnitude 
of the relative Coulomb-nuclear contribution on "ddet is quite smooth. The total contribution of the 
Coulomb-nuclear interference can be neglected with quite a good accuracy. 

4.6 Nuclear contribution to the rotation angle 

The relative contribution of nuclear scattering of deuterons by target nuclei fnn/^o is obtained, using 
the approximate expression (3.7): 



^0 



kRl 



Red 



Tmdi 



Imd) 

Redi J 



1 — exp 



2 + k'^m ^2 



1 



1 



(4.42) 



Similarly to the previous case, let us represent the dependence of (pnn/^Q on the angle of the 
detector for three values of average squared angle of multiple scattering at depth z: 



Vrm. 




Fig. 5: Relative nuclear contribution s as a function of the angle of the detector for 02_|^ = 0.23 • 10 ^ 
cm~^ (solid curve), 6*^2 = 0.23- 10~^ cm~^ (dashed curve) and = 0.23- 10~^ cm~^ (dotted curve). 



In the case of a thick target On, the relation ipnn/vo at large thickness becomes independent 
of i?det ^'iid equals its maximum value. 

f— V = AfRed^-ImdV (4.43) 
The numerical value of \<Pnn/ folmax ~ 0-24. 

For 9z ^ On, the contribution of (fnn/^o depends on i^det (see Fig.5); for -d^et ^ Oz this relation 



is I = — 7T Red— — ; Imd 0'^. It increases gradually with growing angular width of 

\ (fo J kRj V Redi J 2 

the detector and achieves the maximum absolute value (4.43) at -ddet ^ ^n- 

Prom this follows that the nuclear contribution can be neglected for any i?det- In this case 

the magnitude of the relation \'Pnn/'Po\ is of the order of 10~^. 



Conclusion 

Experimental observation of the phenomenon of birefringence of particles naturally always implies 
such experimental arrangement, in which the beam of particles is directed to the target, passes 
through it and then the detector registers the polarization characteristics of particles moving in the 
direction of incidence of the initial beam and detected within a certain small angular interval relative 
to this direction. A comprehensive description of the interaction between particles (dcutcrons) and 
the nuclei of matter is given by a kinetic equation for the density matrix. Using this equation one 
can also see that the scattering process may be both coherent and incoherent. The birefringence 
effect itself, and, in particular, one of its parameters, the rotation angle of the polarization vector 
is the result of a coherent interaction between the particle and the target. Alongside with coherent 
interaction, incoherent interaction also leads to the rotation of the polarization vector relative to 
the chosen direction. That is why the total effective rotation angle is the characteristic, which is 
measured in the experiment. As it has been shown, real experimental conditions also influence the 
value of ipeff- being defined as an integral characteristic, it depends on the value of the angle of 
the detector collimator. In the general case, the magnitude of the Coulomb-nuclear and nuclear- 
nuclear contributions increases with growing i^det- In particular, at 'ddet > On, the correction from the 
Coulomb-nuclear interference exceeds by one order of magnitude the main contribution of ipo. Here 
arises the question about the magnitude of the total contribution of the Coulomb-nuclear interference 
to the rotation angle (feff- It has been shown that these terms are compensated completely only 
when the beam particles are registered in a 47r geometry of the experiment and the deuteron energy 
is such that one can take into account only the first Born approximation to the Coulomb amplitude. 
Therefore, in a real experimental arrangement, the compensation of the Coulomb-nuclear terms will 
not be observed. The quantitative analysis demonstrated that for a thin target, the magnitude of 
the total Coulomb-nuclear interference for 'ddet ^ On equals zero with high accuracy. At the same 
time, as a result of multiple scattering of particles in the target, the contribution of the Coulomb- 
nuclear terms becomes dependent on one more parameter: the average squared angle of multiple 
scattering. It has been shown that if its magnitude at depth z is of the order of or larger than the 
magnitude of the squared diffraction angle of nuclear scattering, the Coulomb-nuclear interference 
can be neglected for any angles of the detector. If this condition is not fulfilled, the contribution of 
the total Coulomb-nuclear interference to the rotation angle should be taken into account. 

The contribution of the nuclear-nuclear interaction for the two cases of target thicknesses consid- 
ered above is small as compared to the main contribution. 

Thus, the elastic coherent scattering of deuterons in the target leads to additional corrections 
to the birefringence effect, the magnitude of these contributions appreciably depends on both the 
average squared angle of multiple scattering and on the specific geometry of the experiment, namely 
on the angular width of the collimator. 
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